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Abstract. The European power grid can be divided into several market areas where the price of electricity 
is determined in a day-ahead auction. Market participants can provide continuous hourly bid curves and 
combinatorial bids with associated quantities given the prices. The goal of our auction is to maximize the 
economic surplus of all participants subject to transmission constraints and the existence of linear prices. 
In general strict linear prices do not exist in non-convex markets. Therefore we enforce the existence of 
linear prices where no one incurs a loss and only combinatorial bids might see a not realized gain. The 
resulting optimization problem is an MPEC that can not be solved efficiently by a standard solver. We 
present an exact algorithm and a fast heuristic for this type of problem. Both algorithms decompose the 
MPEC into a master MIP and price subproblems (LPs). The modeling technique and the algorithms are 
^\ applicable to all MIP based combinatorial auctions. 
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1. Introduction 



Energy markets (and in particular electricity markets) in Europe are required to be liberalized due 
to a directive of the European Union. While this liberalization is a step toward fostering competition, 
efficient prices, and reduced volatility, new mechanisms have to be put into place to realize this process. 
In particular markets have to be defined, implemented, and run in order to allow for the trading 
of energy. In order to guarantee a smoothly running, efficient market the correct market design is 
imperative. Whereas markets for investment commodities (and investment assets in general) often have 
a rather simple structure, markets for consumption commodities have to account for the particularities 
of the traded commodity, in this case electrical power. For example, electrical power can only be 
stored to a very small amount and energy resulting from overproduction has to be properly disposed 
of. Both procedures are rather expensive so that under certain circumstances it can be economically 
feasible that electrical power can have a negative price, i.e., the buyer is paid for buying the electricity. 
Moreover ramping up and shutting down a power generation facility can take a significant amount of 
^—i time (several hours up to days) and so energy providers might be discouraged to trade larger amounts 

of energy due to overproduction resulting from extensive shutdown periods. To counterfeit this trend 
and to ensure a sufficient market liquidity which is imperative, many electricity markets allow for bids 
which exhibit a combinatorial structure. For example, an energy provider might be willing to offer 
a certain amount of electricity in the morning if and only if somebody will also buy the energy that 
is produced (inevitably due to long shutdown phases) in the afternoon. In the absence of such bids, 
the energy provider might not have any incentive to offer electricity in the morning as the profit in 
the morning might be eaten up by the disposal costs in the afternoon resulting in a net loss. These 
combined bids are executed en bloc. Electricity exchanges (e.g., APX, Belpex, EPEX Spot, Nord Pool 
Spot AS) implementing these structures have been successfully established since the early 90s. 

Markets with such an enhanced structure are commonly referred to as smart markets and clearing is 
performed by solving the underlying (often combinatorial) optimization problem. In contrast to more 
mundane market designs that are cleared via order book uncrossing, typically smart markets require a 
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significantly larger amount of computational resources and methods from combinatorial optimization. 
The electricity market that we consider here is such a smart market given e.g., its block bids and flex 
bids. Market participants can bid for electricity by submitting bid curves, block bids, and flex bids to 
their national exchange. These bids specify (for every hour) the quantity (or combinations of those) 
that the participant is willing to buy or sell for a given market price. 

In a consequent next step the regulators additionally suggested a coupling of these liberalized na- 
tional markets with each other. This international cross-border trading between different electricity 
markets allows for leveling prices in adjacent market areas leading to overall more efficient markets 
and removing arbitrage opportunities that could be exploited by energy companies operating in differ- 
ent market areas at the same time. The current regulation on cross-border electricity exchange (cf. |EC| 
1228/2003 [2003]) of the European Union enacted in July 2003 requires a market based allocation of 



transmission capacity of congested cross-border transmission lines. A study (cf. Consentec and Fron 



tier Economics [2004]) commissioned by the European Commission concludes that the allocation of 
transfer capacity should be done via an implicit auction. 

In this article we present the market design that is currently in use by most of the European electricity 
exchanges. Basically it is a welfare maximizing combinatorial auction that determines linear clearing 
prices. Even though the main focus of the paper will be on solving the underlying problem, we also 
explain why the market model evolved to its current state. Furthermore we will discuss the differences 
to unit commitment models with nonlinear clearing prices. 



The US Market At first we want to give a short introduction to electricity market clearing models that 
are used in the US. In the US there exist many Regional Transmission Organization Markets (RTO Mar- 
kets). Theses markets are run by Independent System Operators (ISO) who are using unit commitment 
(UC) models to determine a cost minimal dispatch. A standard reference for unit commitment models is 



Hobbs et al. [2001 ]. These models allow for submitting the cost structure and production capacities of 



each power plant. In particular it is possible to submit startup costs and minimal production capacities 
which are modeled by using binary variables. This detailed information allows the ISOs to determine a 
cost minimal commitment of each power plant. Under the assumption that all constraints that describe 
the unit commitment problem are linear and that no binary variables exist. Then the optimal dispatch 
can be determined by solving an LP and the prices that support this solution are given by the dual vari- 
ables. Such prices are called linear prices, because the function that prices the commodity a depending 
on the quantity q is a linear function (_f a (q) := n a q). As soon as binary variables are involved we have 
to solve a mixed integer linear program (MIP) and the determination of reasonable prices that support 



an optimal solution is not as easy as in the LP case. O'Neill et al. [2005 ] present a pricing scheme that 



allows us to define prices for a MIP solution. The main idea is, that we have to add the binary actions to 
the commodity space. Then we can find linear prices for all commodities, including the binary actions. 
There already exist auctions that are using similar pricing schemes, for example the New York Indepen- 
dent System Operator (NYISO) and the Pennsylvania-New Jersey-Maryland Interconnection (PJM). Even 
though we could use a similar approach for the European market the participants would not accept the 
prices as they are discriminatory: every participant who submits a combinatorial bid would receive or 
pay an individual side payment and the auctioneer would keep the entire economic surplus generated 
by the combinatorial bids (cf., Appendix^). Another technique that allows us to define prices for a MIP 
solution is called convex hull pricing (Gribik et al. [2007], Wang et al. [2011 ]). This pricing scheme 
reduces the impact of side payments but still requires the use of individual side payments yielding 
discriminatory prices. 



The European Market. Similar to the US there exist several regional Transmission System Operatorqj 
(TSO) in Europe. The main difference to the ISOs in the US is, that the European TSOs only maintain 
system feasibility and reliability. They do not control the spot market for electricity, but they determine 
network boundary conditions within which they can guarantee the system feasibility and reliability. So 
the boundary conditions (e.g., available transmission capacities) are submitted to independent power 
exchange^ who will run the financial spot market auctions subject to the network constraints. On top 
of these power exchanges a market coupling system^] collects all order book data from the exchanges 
and network constraints from the TSOs and computes economically efficient power flows between the 
power exchanges. These flows and price signals are submitted back to the exchanges and TSOs. Finally 
the exchanges run their own auction incorporating the cross-border flow signals and the price signals. 
All power exchanges that we consider here determine linear nondiscriminatory clearing prices for 
their specific bidding area. This implies that the market coupling system must also determine linear 
nondiscriminatory clearing prices, and thus we can not use a pricing scheme that is similar to the one 



proposed in O'Neill et al. [ 2005 ] . We have to satisfy the pricing rules that are given by the exchanges. 
As of now, the TSOs do not submit flow based network constraints (PTDF-ma trices) . They only submit 
upper and lower bounds for the flows and ramping conditions. This is due to the fact that the coupling 
algorithms used within the Nordpool Spot region and the CASC-CWE region do not use PTDFs yet. 
Therefore a market coupling system that couples these two regions cannot incorporate PTDFs, however 
the algorithms that we present can be easily adjusted to solve a flow based market coupling problem. 



The electricity market design that we consider here is based on Kirschen and Strbac [2004] and 



Hogan [1998]. It uses bid curves representing the marginal cost of production and the marginal ben- 
efit of consumption and prices the congestion of network branches. It is known that under certain 



circumstances this market design might be exploited by strategic bidders (cf. Woo et al. [2003]) and 



there are even indications that this might have been the case at the German electricity exchange EEX 
(cf. Canty [2009]). Whereas this might be a real-world issue that has to be addressed by regulators, 
we assume a smoothly functioning market as our focus is a different one. There exist several articles 
addressing the particular issues of clearing energy markets however with different models: Most of the 
models either do not contain binary decision variables or they use a nonlinear pricing scheme. For ex- 



ample, the model introduced in Kirschen and Strbac [2004] considers only a single time slot and only 
continuous hourly bid curves. While this fundamental approach provides a very nice illustration of the 
equilibrium conditions and security constraints it does not provide a solution approach in the presence 
of combinatorial structures. The approach in Haubrich et al. [ 2006 ] is a market coupling model with 
step functions as bid curves and without the fill-or-kill condition of block bids, i.e., block bids can be 
partially executed in this model. The resulting model is a linear program that can be solved by any 
standard solver. The most difficult part however, the ensuring of the fill-or-kill condition in combination 
with linear prices is not addressed by this approach. On the one hand introducing the fill-or-kill condi- 
tion will include knapsack problems leading to NP-hard optimization problems. On the other hand the 
inclusion of binary variables enforces us to manually ensure the existence of linear prices. This is done 
by adding additional optimality conditions to the constraints. The resulting model will be a Mathemat- 
ical Program with Equilibrium Constraints (MPEC) . It is very difficult to handle these complementarity 
conditions in conjunction with the binary decision variables. The preliminary approach presented in 
Krion [J2008] handled the fill-or-kill conditions as well as the existence of linear prices with respect to 
the three basic bid types. Therefore the bid curves where modeled by the delta method with the help of 
additional binary variables. However the complementarity conditions of cross-border-trades where not 
included in this preliminary model such that the fill-or-kill conditions could cause adverse flows from 



TSOs within the considered market: Amprion, CASC-CWE, elia, EnBW, Energinet.dk, Fingrid, 50hertz transmission, Rte, 
Statnett, Svenska Kraftnat, Tennet, Transpower, ... further European TSOs can be found at www . entsoe . eu 



2 Power exchanges within the considered market: APX, Belpex, EPEX Spot, Nordpool Spot 
3 In our case the system is run by European Market Coupling Company (EMCC) 
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a high price area to a low price area. The stakeholders decided that adverse flows may only be caused 



by transmission constraints. This requirement was added to the model by using a heuristic (see [ EMCC 
GmbH| Chapter 4.3.2] and Section[7]for details). 



Our contribution. In this article we use the market model currently used by European electricity 
exchanges. The presented model is compliant with the European regulatory framework and approved 
by the regulator. The European market model differs from common US RTO market models that are 
using nonlinear pricing schemes. The model incorporates binary decision variables and guarantees the 
existence of nondiscriminatory linear clearing prices. The existence of linear clearing prices is enforced 
by complementarity conditions yielding a Mathematical Program with Equilibrium Constraints (MPEC). 
We present a deterministic algorithm yielding a provably optimal solution and a fast heuristic that 
exploits the model structure. Instead of solving the complicated MPEC directly we decompose it into a 
Mixed Integer Quadratic Program (MIQP) and a linear pricing problem. Empirical tests suggest that the 
solutions determined by the heuristic are optimal in many cases. We also ensure unique prices without 
changing the economic surplus. Additionally, we introduce rules that allow market participants to check 
whether the determined prices satisfy necessary optimality conditions. 

Outline. The structure of the article is as follows. In Section [2] we give a brief introduction to the 
European electricity market and its features. In Section [3] we introduce the combinatorial optimization 
problem that is used to determine optimal flows and efficient prices. Analyzing this model we obtain 
optimality conditions which are formulated in Section [4J We then establish in Section [5]the uniqueness 
of prices and then we formulate a clearing heuristic and an optimal algorithm in Section[6J We conclude 
with some computational results in Section [7] and some final remarks in Section [8] 

2. Preliminaries 

In the following section we recall the structure of the electricity markets and its features. We also in- 
troduce the necessary notions and notations. We first explain the electricity market related terminology 
and then, in the second part, we will introduce the notation that we will use later. 

European market design. The design of the day-ahead electricity market that we consider here is 



an extension of the designs of Kirschen and Strbac [ 2004 ] and Haubrich et al. [ 2006 1 and meets 



the needs of Central West European Market Coupling Project (CWE MC) and Nord Pool Spot AS (NPS). 
Nonetheless the model can be easily adapted to account for particularities of other electricity exchanges, 
if the regulatory boundary conditions are similar to the European regulatory conditions. 

In a day-ahead energy market participants bid today for electricity that is being bought or sold on the 
following day. In fact, participants bid for electricity in specific time slots on the next day. In a coupled 
day-ahead electricity market, the participants can bid for electricity in and across certain market areas. 
The participants are allowed to submit their bids up to a predetermined point in time. We basically 
differentiate between three types of bids that will be explained in more detail later. The first type is 
an hourly bid curve where the participant determines the amount of energy that he is willing to buy or 
sell given a certain price. The second type of bid are so called block bids. Using this type of bid the 
participant can bid for block structures of energy, e.g., selling a certain amount of energy in the morning 
and a certain amount of energy in the evening. These bids have a fill-or-kill character as they can only 
be executed as a whole and exactly this type of bid introduces combinatorial features into the market. 
The last type of bids are so called flexible bids (short: flex bids). These bids allow for selling a certain 
amount of electricity in exactly one time slot. The actual slot in which this takes place is not specified. 
All bids are valid for specific time slots (except for flex bids) and specific market areas. As soon as 
the submission deadline has passed the auction algorithm has to determine an optimal matching of 
buyers and sellers, determine efficient prices for each market area and each time slot, and compute the 
optimal cross-border flows between adjacent market areas. The determined price reflects the price for 
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Figure 2.1. Modeling of the hourly net curves and the economic surplus. 



one electrical power unit that is delivered throughout the specified time slot. Although these prices are 
only valid for those bids that were cleared by the auction algorithm, they do also serve the purpose of 
being price signals for the intraday electricity market. 

In order to refer to a specific area we write a s A where A denotes the set of areas and similarly 
we write teTto refer to a specific time slot where T denotes the set of time slots; here hours of the 
day. The price in area aeAat time t e T will be denoted by n at e [p,p] where [p,p] denotes the 
feasible price interval with p <p. We now define the types of bids in a more formal way that suits the 
later exposition. In the following we will consider aggregated demand and supply curves without loss 
of generality. 

Hourly demand and supply curve. Participants can submit an hourly bid curve for every hour and area 
to specify the feasible quantities that will be bought or sold depending on the price. These curves are 
piecewise linear, monotonically decreasing, and are defined on [p,p]. Positive quantities correspond 
to electricity that is bought. To determine the surplus and the prices it is sufficient to regard only one 
aggregated hourly sell and one hourly buy curve per area and hour as we do not discriminate between 
different buyers or sellers. We can further aggregate these buy and sell curves to one combined hourly 
net demand curve per area and hour (see e.g., Figure |2.1|). The net curves contain all necessary 



information to maximize the economic surplus and to determine efficient prices (see e.g., Meeus et al. 



[2005]) on the one hand and facilitates the actual market clearing on the other hand. It is easy to see 



that the hourly net curves are also piecewise linear and monotonically decreasing. We refer to H c N 
as the set of segments. Each segment h e H at C H is attributed to a market area a and a time slot t 
and its execution state is modeled by the continuous variable 5 h e [0, 1]. 

Block bids. Block bids enable traders to buy or sell electricity for several (not necessarily consecutive) 
hour slots with a single bid. We denote the set of block bids by B. Block bids can only be executed 
as a whole (fill-or-kill condition) and the price determined for executing a block bid beB will be 
the weighted average of the prices of electricity in the corresponding areas weighted by the traded 
quantities. Note that even if the determined prices seem to allow the execution of a block bid it might 
not be executed as the execution could imply a loss of economic surplus or the execution could change 
the prices, so that this block or other blocks cannot be executed. These blocks are called paradoxically 
rejected blocks. Whether or not a block bid is executed will be denoted by the binary variables /3 b with 
b e B. As we will see later, finding the optimal block bid selection is the most difficult task faced 
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Figure 2.2. Example of an interconnected electricity market (Source file from Wikipedia). 



when clearing the market. The trader further has the possibility to link certain block bids, so that the 
execution of some block b 1 e B implies the execution of some other block b 2 ^B. 

Flex bids. With flex bids a trader can buy or sell electricity in exactly on hour without specifying the 
hour in which the bid has to be executed. A buy bid can be executed in an hour where the clearing 
price is lower than or equal to the bid price (and similarly for selling electricity). Flexible bids must 
satisfy the fill-or-kill condition as well. We refer to the set of flex bids by F and for each / € F we have 
binary variables ipf >t with t ^T. These variables will be used to determine when a flex bid is executed. 

Cross-border trades. Trading electricity across areas is called a cross-border trade. The cross-border 
trader acts as a consumer in one and as supplier in another area. Cross-border trades are only possible 
if there are interconnectors connecting the areas. The latter will be denoted by C and the transmitted 
electricity will be denoted by t c t with c ^C and t e T. These interconnectors have limited transmis- 
sion capacities so that capacity cannot be allocated arbitrarily. In an implicit capacity allocation, there is 
only one cross-border-trader per interconnector and the transmission quantity will be determined by a 
market coupling system with aim to maximize the economic surplus. In this setting traders cannot di- 
rectly trade across market segments: the cross-border trader at the involved interconnector determines 
the execution quantities by maximizing the total economic surplus. Note that also the surplus of the 
cross-border trader, the so called congestion rent, is part of the economic surplus. 

Notation and definitions. For the ease of presentation we now introduce and summarize the notation 
that we will use in the remainder of the article. 
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Sets. The following sets are used in the model: 

T c N hours of the day, ordered 

A market areas 

H c N nodes of all hourly net demand curves 

H a >t c H nodes of hourly net demand curve in area a, hour t 

S c H curtailment segments of hourly net demand curves 

B block bids 

B a c B block bids in area a 

F flexible bids 

F a d F flexible bids in area a 

C c A x A interconnectors 

CJ" = {(r,s) e C | r = a} interconnectors starting in a 

C~ = {(r,s) G C \s = a} interconnectors ending in a 

IcBxJJ links between block bids 

Variables. We use the following variables in the optimization problem: 

clearing price in area a e A in hour t sT 

(1 block bid b e B is executed 
n 
no execution 

J 1 flexible bid / e F is executed in hour t e T 

1 no execution 
partial execution of the hourly net demand curve segment 
with the index heH 
T c,t G [^ c t>^c,t~\ fl° w on interconnector c = (r,5) e C at time teT 

T ct > denotes a flow from r to 5 at time t 

Parameters. We now introduce the parameters that we use in the optimization model. The parameters 
used in the optimization model are dimensionless. In the actual implementation quantities are given in 
megawatt (MW) and prices are denoted in Euro/MW A positive quantity denotes a demand whereas 
a negative one indicates supply. All submitted price limits are required to be in the price interval 



n a,t 


e [p,p] 


Pb 


£{0,1} 


Vf,t 


€{0,1} 


5 h 


6 [0,1] 



P = [p,p] [Nord Pool Spot AS 2011 Section 2.3]. This implies that we can find prices within this 
interval, thus the constraint n at e P is not a restriction. 

Parameters of block bids. Let b e B a be a block bid in area a e A. Then 

q b >t e R denotes the quantity in hour tel and 

p b e P denotes the average limit price. 

The average limit price p fe is the average price, weighted by the traded quantities. More precisely, the 
block bid b can only be executed if the trader can realize a net profit, i.e., Xlter^Pfe — n a,t) c lb,t — 0- 

Parameters for linking block bids. Let b SB and ceBbe block bids. The set L C B x B denotes the links 
between block bids. If (b, c) e L, then b can only be executed if also c is executed. 
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Parameters of flexible bids. Let / e F a be a flexible bid in area a e A. Then 

q_f e R denotes the quantity and 

Pf e P denotes the limit price. 

The flexible bid / can be executed in exactly one hour t e T if the trader can realize a profit when 
executing in t, that is if {pf — n at )qf > 0. 

Parameters of interconnectors. Let c = (r,5) e C be an interconnector and let t ^T. The interconnector 
is a directed edge, whereas r denotes the source and s denotes the sink. Nonetheless, the actual 
transfer on this edge is bidirectional and the transfer direction is signified by the sign of t c t . We use 
the convention that a negative z ct indicates transmitted electricity against the direction of the edge, 
i.e., from s to r. Furthermore, the volume transmitted via an interconnector might not be allowed to 
change too much between two periods. The allowed absolute difference (between consecutive slots) 
is the ramp rate. There are also upper and lower bounds on the transmission quantity called available 
transfer capacities (ATC). We write 

T y €l as the upper bound for the flow on c (ATC), 

t c t e E as the lower bound for the flow on c (ATC), 

T c _i el as flow on c in the last time slot of the previous day, 

t c G Rt U {00} as the ramp rate. 

We put t c = 00 if we have no ramping on c. 

Parameters of hourly net demand curves. We consider one aggregated hourly net demand curve per area 
a and time t. An hourly net curve is a piecewise linear function in M 2 consisting of (price, quantity) 
tuples. We require that the function is monotonically decreasing, i.e, whenever we have two feasible 
tuples (pi, qi), (p 2 , Q2) with p 1 < p 2 then q 1 > q 2 holds. This condition ensures that the net curve 
is well defined, i.e., for an increasing price the net demand decreases (not necessarily strictly). Every 
point of the curve represents a feasible (price, quantity) combination and as a whole they constitute the 
feasible price range depending on the traded quantity. 

Let h e H at be a defining node of an hourly net curve in area a at time t. Then 

q e M~ denotes the minimum net demand quantity, 

Pf,eP denotes the price coordinate of the node, 

qi,€l denotes the quantity coordinate of the node, 

Ap h €E Mg denotes the price difference to the next node, 

Aq h €E Mg denotes the quantity difference to the next node. 

The hourly net curve must intersect the price axis and has to be defined on the entire feasible price 
interval P. If the net curve has no intersection with the price axis, then a vertical segment is added, 
that connects the price axis with the curvejj This segment models the curtailment of hourly bids and 
is added to the list of curtailment segments S. Effectively, it sets a maximum price for which electricity 
might be bought. An example of a net curve arising from a demand and a supply curve are depicted 



in Figures 2.1 and 3.1 It is possible that quantity or price do not change over consecutive intervals. 
Without loss of generality we assume that the function remains constant only over exactly one interval 
at a time, i.e., there are no three tuples (pi,qi), (P2> ( 3 , 2)> an d (£3.33) with either pi < p 2 < P3 and 
Qi = Q2 = Q3 or <h — ?2 — ^3 an d Pi = Pi = P3; i n this case we could simply drop {jp^q?). This 

Note that this is a curtailment rule of CWE MC, that is used in cases of non-harmo nized price rang es. In cases of 



harmonized price ranges curtailment is shared among adjacent areas. For more details see [EPEX Spot 2010 Chapter 6.5.2] 



assumption simplifies some proofs on the one hand and, when used within the preprocessing, speeds 
up the solution process on the other hand. 

3. Model 

In this section we will present the complete MPEC that models our clearing problem. We first intro- 
duce the overall model structure and then discuss each constraint as well as the objective function in 
detail. 

max (economic surplus) J] \(p h + Ap h ) Aq h 5 h - -Ap h Aq h 5^ + J]p b q b ,t Pb + ^Pf If <Pf,t 

heH beB feF 

t£T teT 

(MPEC) 

s.t. (quantity constraints) 

[n] Va e A, t e T i a t + X Aqh5h + X qh ^ h + X q f ^U = 2 T ^ ~ X Tc > f 

heH a>t b£B a f£F a cec~ ceC+ 



(Cll 



[p]VceC, t€T - t c < T C)t - icj-t < t c dC2| ) 

V(b,c)£L /3 b <l3 c (Q 

VbeB /SbG{0,l} 

V/eF,ter V/>t e{0,l} 

[pi] Vc e C, t e T T ct e [T ct , r ct ] 

[v] VfreH 5 h e [0,1] 

(price constraints) 



VaeA,/eF a ,ter ^(py - n a>t ) sgn(q f ) > flC6D 

VceC,teT (M Cj t - M c t ) + (Pc,t - Pc,t) + [(Pc,t+i - Pc,t+i)] + re r,t - ^s.t = flC7l ) 

VceC,teT (T C)t - T ct ) p ct = flcTI ) 

VceC,teT (-^c,t + l c ,t)M ct = flC7l ) 

VceC, reT (T ct - T C;t _! - t c ) p ct = flC 7\ 

VceC,ter ( T c,t-i _T c,t _;? c)Pc,t = dC 7} 



VceC,ter jLt ct , ju^ t , p c<t , p Cit > ( |C 7D 



-c,t 

•* 



VaeA,ter,heH a]t - Aq hPh (5*) + 7i a]t Aq h + v h - v h = ((Cj) 

VaeA,ter,heH at (5 h -l)v h = JC8] ) 

VaeA,ter,l!eH a]t (-<5,Jv h = flC~8] ) 

VaeA, teT,heH at v h , v h > JC8] ) 

VaeA, t€T 7i at e P 

Note that the price conditions contain complementarity conditions. The variables p,/i, and v are dual 
variables of the individual optimization problems of price takers. As the complete model incorporates 
the individual optimization problems, these variables are primal variables in the MPEC. We need to 
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explicitly add the price conditions to the model, because they can not be fulfilled in general if they are 
not included in the model. 

The main idea of the model is to define prices that are consistent with the dual prices that we know 
from continuous models: Assume that we solve a model, that only contains the quantity constraints 
(except for constraint < \C 3] )) and that the model has the integral property. Then we can find dual prices 
to an optimal primal solution that satisfy the price conditions. In this case, the price conditions could be 
added to the model, without cutting of any optimal solution. This assertion can be verified by writing 
down the KKT conditions of the QP-relaxation. 

3.1. Constraints. We now introduce the constraints mentioned above and provide justification and 
interpretation of those. 

C 1. clearing condition. This constraints ensures the identity of executed net demand and net 
import. The constraint is similar to the classical conservation of flow constraint typically present 
in network flow models. In every area and every hour the executed net demand must coincide 
with the net import: 

VaeA, teT Q a t + X A ^ 5 ft+ J] <lb,tPb + ^<lfVf,t= X T ^~ X Tc > f • 

heH at beB a feF a ceC- ceC+ 

C 2. ramping condition. This ramping condition incorporates the ramping requirements on the 
interconnectors by limiting the change between two consecutive time intervals accordingly. 
From one hour to the next hour the flow on interconnector c may only change by t c units: 

VceC, teT |T Cjt - T Cjt _:iJ < t c . 

C 3. block bid links. To model the linking of block bids we add the following constraint. If (b, c) e 
L, then block b may only be executed if block c is executed: 

V(b,c)£L /3 b <l3 c . 

C 4. flexible bid execution at most once. A flex bid / can only be executed once in one specific 
time slot: 

V/eF X! 1 ^- 1 - 

C 5. block price condition. This constraint ensures that our block selection strategy is incentive 
compatible. In principle it is possible to obtain an equally high or even higher economic surplus 
by requiring altruistic behavior from some of the participants. However this is not in line with 
general economic principles as no participant is willing to sell at a loss. Therefore, a block bid 
can only be executed if it does not generate a net loss: 

VaeA,beB a p b ^Ob - ^a,t) <lb,t ^ ° • 

teT 

This is a quadratic non-convex function that is difficult to handle, but it is easy to linearize it if 
we introduce a big constant M b for each feeB. We know that /3 b is a binary variable and we 
know that p e [p,p]. Now we can rewrite the constraint as follows: 

Y^(Pb- n a,^qb,t >M b (l-/3 b ) with M b :=X min KPb-70qb, t I ^ e [poli- 
ter teT 

C 6. flexible bid price condition. Moreover, the execution of a flex bid is contingent on economic 
compatibility. More precisely, a flex bid is only executed if it does not generate a loss: 

Va e A,f e F a , t e T Vf, t (Pf ~ n a,t) sgn^) > 0, 
where sgn(.) denotes the sign function. 
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ceC teT 




[p] Vc e c, t e r 


T c ,t ^ T C)t 


[ju] Vc e C, t e T 


- T ct < "let 


[p] Vc e C, t e T 


T Cjt -T C;t _i <T 


[p] VceC,teT 


T c ,t-1 - T, t < T 



C 7. flow price condition. This condition represents the optimality condition of the cross-border 
traders optimization problem. The regulators force them to be price takers such that the prices 
are externally given and the only remaining decision variable is the flow T ct . The individual 
optimization problem is to maximize the congestion rent subject to the available transmission 
capacity and the ramp rate: 

max ^^(Jt, lt -tr,t)T Cl t (LP-TSO) 

ceC teT 
s.t. [jS] Vc S C, t e T T r t < T r t ATC 

ATC 

ramp rate 

ramp rate 

Basically an optimal flow implies that there may only be a price difference between two ad- 
jacent areas if at least one transmission constraint is active. In other words, the prices of 
adjacent areas coincide, if there is free transmission capacity between those areas. From now 
on the KKT-conditions of ( |LP-TSOD will be called the flow price condition. 

Definition 3.1 (flow price condition). Let t S R CxT , t c t < t < t c t , and n S P AxT . The flow- 
price tuple (t, tt) satisfies the flow price condition if there exist ju, ju, p, p e l CxT , so that for all 
c = (r,s) eC and t eT 

Mct.M^.Pc.t.Pct >0, (3.1) 

(ju ct -/*.) + (Pc,t ~ Pc,t) + [(Pc,t+i - Pc,t+i)] + Tlr.t - ^s,t = > C3.2) 

( T c,t - T C;t ) ju c ,t = ° - C3.3) 

(-T c , t +l C;t )M ct = 0, (3.4) 

( T c,t-' r c,t-i- ;? c)Pc,t = , C 3 - 5 ) 

(T c ,t-i ~ T Cit - T c ) p Cit = . (3.6) 



7n terms of duality theory the conditions fl3.1[ ) and ( |3.2| ) correspond to dual feasibility conditions 
and the remaining ones are complementarity conditions. The term in square brackets is only added 
ift + leT. 

If there is no ramping on an interconnector c = (r,5), this condition simplifies to[j 

(< , T r f < T r t 
~ c ' 1 c ' f . (3.7) 

>0 ,T C;t <T C;t 

The economic interpretation of this simplified condition is that free transfer capacity in the 
forward direction implies that the price n rt in the source area must be higher or equal to the 
price K st in the sink area. If this would not be the case, then the free capacity should be used 
to transfer electricity from the source to the sink, and therefore the market would not be in an 
equilibrium. 

If, on the other hand, the ramping condition (C 2) on an interconnector c = (r,s) is only 
binding (i.e., satisfied with equality) in two consecutive hours t and t + 1, and there is still 
available transfer capacity in that time window, then the condition simplifies to: 

n s,t - K r,t = -(X,t+1 ~ K r,t+l) ■ ( 3 - 8 ) 



A derivation of the equations l |3.7| l to ( |3.1Q| > can be found in Appendix B 
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This establishes that the price difference at one point in time can have an impact on the price 
difference in the following ones. The condition also implies, that prices of adjacent areas 
can only deviate if at least one respective transmission constraint is active. In particular it is 
sufficient, that in one hour the transmission is limited only by the ramping condition and not 
by the available transmission capacities to trigger a price difference. It holds (terms in square 
brackets are only present if t + 1 e T): 



Vc = (r,s)<EC, te T 



n r,t > n s,t 



L c,t 
*c,t = T 



T Cjt V 



T r t -T 



c,t 



V 



■c,t-l 



C,t— 1 — l C 

T r t T r 



[v 
[v 



T r t -T 



c,t+\ 



•c,t+l 



T t =T c ] . 



(3.9) 
(3.10) 



C 8. filling condition of hourly net curves. This constraint represents the optimality conditions 
of an hourly bid. Participants who submit hourly bids are modeled as price takers. In the 
optimization problem associated with an hourly bid h the price is exogenously given and the 
only remaining decision variable is 5 h . The objective is to maximize the economic surplus. In 
the next chapter we will explain the following integral in more detail. 

r 5 h 
max Aq h (p h (u) — n at )du+ const (QP-Hourly) 



s.t. [v] +5 h <l 
[v] - 5 h < 

with P/i(<5/j) = p/j + (l — S/jAp/j. From now on the KKT conditions of this optimization problem 
will be called filling condition. 

Definition 3.2 (filling condition). Let 5 S [0, 1] H and n S P AxT . The tuple (5, n) satisfies the 
filling condition if and only if there exist variables v, v e M. H , so that for all a e A, t e T, he H at 

- Aq h p h (5*) + n a>t Aq h + v h -v h = 

(5 h -l)v h = 

(-5 h )v h = 

This condition can be reformulated to an equivalent condition: 

Definition 3.3 (filling condition). Let 5 S [0, 1] H and n S P AxT . The tuple (5, n) satisfies the 
filling condition if and only if 



Va <E A, t e T, h e H a t with Aq h ^ 



with p h (5 h ) = p h + (1 - 5 h )A Ph . 



7i„, i =p h (5 h ), if 5 h e (0,1) 
>p h (5 h ), if5 h = 



(3.11) 



(3.12) 



The second definition directly shows that an optimal solution only contains quantity-price- 
combinations (5, 7i) that were defined by the piecewise linear bid curves. A MIP formulation of 
the filling condition based on the delta-method (cf. Dantzig [ 1998[ ]) can be found in Appendix[c] 
and a graphical illustration is given in Figure |2.1 



For the sake of completeness we want to list two further constraints that are relevant in extreme situa- 



tions but will not be handled by the model ( ]MPEC| ) . We omitted these constraints so that we can focus 
on the constraints listed beforehand. 
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Figure 3.1. Extending hourly net curve from small price interval P a to P. 



Remark 1. Combinatorial bid curtailment. This constraint is a market rule requested by the ex- 
changes to improve the system stability: In curtailment situations price indifferent hourly bids have 
a higher priority than combinatorial bids. If the hourly net demand is curtailed in area a at time t, 
then combinatorial bids with demand in area a, time t will also be curtailed (i.e., they must not be 
executed). Similarly this should hold for curtailed supply. The set of net curve segments that represent 
the curtailment is denoted by S. We have the following condition: 



Va<EA,t<E T,h<EH at r\S 
Demand curtailment: 

Supply curtailment: 



If q h > 0, then 
If q h < 0, then 



'Vb eB„ withq b t >0 
V/ e F a with q f > 

r \/b eB a withq b t <0 
yf e F a with q f < 



Pb<5 h 
<Pf,t ^ 5 h ' 

Pb < a - 5 h ) 

ip f ,t < (1 - 5 h ) 



Remark 2. Different price intervals. The power exchanges that need to be coupled have different 
feasible price intervals. The prices in the Nordpool Spot region are limited to [—200,2000] €/MW 
whereas the prices in the CASC-CWE region are limited to [—3000,3000] €/MW. Suppose that the 
feasible price interval P a of area a is a proper subset of P. In this case the hourly net curve is a priori 
only defined on P a . We therefore need to extend the curve to cover the whole price interval P and to 
ensure the existence of prices with the above mentioned properties in P. This is done by extending the 
curve in a horizontal manner. If the curve has no intersection with the price axis, i.e., the maximum 
supply does not cover the minimum demand or vice versa, then the curve is extended in a vertical 
manner at first to create an intersection with the price axis. This vertical segment is called curtailment 



segment and will be stored in the set S c H. An example is shown in Figure 3.1 where the extension is 
indicated by the dotted lines. 

If a solution contains a price n at ^ P a , then n at is changed to the nearest price in P a . By doing so, 
we might violate the flow price condition. In other words, there might be a price difference between 
two adjacent areas although none of the transport restrictions is binding. In fact, this might happen if 
one of the areas has a small price interval and reaches the maximal possible hourly demand or supply. 

3.2. Objective Function. The economic surplus is the sum of consumer surplus and producer surplus 



(cf. Pindyck and Rubinfeld [2005]). It is also called global welfare (cf. Kirschen a nd Strbac [2004]) or 
social welfare. For example a consumer buys q units at a price of n per unit and he is willing to pay 
a price of p > n, then his consumer surplus is given by (p — n)q. The producer surplus generated by 
hourly bid curves is determined accordingly by integrating the clearing price minus the producers cost 
over the quantity from zero to the produced quantity. As we are only interested in the maximization of 
the economic surplus, it is sufficient to compute the economic surplus of the hourly net demand curve 
and the other market products instead of using the hourly demand and supply curves. 

We will now define the surplus of each market product. These definitions are consistent with the 
above mentioned one. Beforehand we define the incremental price and quantity functions in a more 
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Figure 3.2. Economic surplus of an hourly net curve segment. 



precise way. These functions correspond to the linear interpolation between two consecutive defining 
nodes. 

Definition 3.4. For area a e A, hour t e T, and segment h € H at , the net curve segment functions 
p h , q h : [0, 1] -> E are given by: 

p h (z):=p h + <il-z)Ap h , 

q?t(z) := % + Aq h z , 



with q 



min{0,q h -AqJ ifq h >0 
—Aq h , otherwise 



Note, that if Aq h ^ 0, then the inverse supply /demand function is given by Ph(q ft 1 (u)) where u e 
[q h ,q + Aq h ] is the net quantity that is contributed by the segment h. With this definition we can 
define the surplus for net curve segments and hence for the whole curve by additivity 

Definition 3.5. The surplus of a non-horizontal net curve segment h e H at in area a and hour t with 
Aq h j^O is given by 



co h (5, n) := (Phi.% V)) ~ n a,t) du 

Jo 



5 h 
■■Aq h I p h (z)dz - n a>t q h (5 h ) - Aq h 
/o 



rq^(o) 



PhO)dz 



(3.13) 



(3.14) 



Note that the last term in Equation 3.14 is a constant. Horizontal segments do not generate a surplus, as 



they only represent a synthetic bid with quantity zero: For h e H with Aq ft = 0, we have co h (5, n) = 0. 



Expression Q3.141 ) is equivalent to the net consumers/producers surplus in [Kirschen and Strbac 2004 
Chapter 2.2]. 

Note that the surplus for the net curve corresponds to the area between the curve and the price 



axis. Examples are provided in Figures 2.1 and 3.2 Next we define the surplus for the combinatorial 
structures, i.e., the block and flex bids. 

Definition 3.6. The surplus of a block bid b e B a with a &Afor n e P AxT and /3 G {0, 1} B is given by 



co b (n, /3) := Y^Pb ~ n a ,t)<lb,t Pb- 



t£T 



Definition 3.7. The surplus of a flexible bid /eF a with a ^Afor n e R AxT and (p e {0, l} Fxr is given 

by 

co f (n, y) := Y^Pf ~ n a,t)qf <P fit- 



ter 
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Another surplus arises on the interconnectors. This surplus represents the surplus of cross-border 
trades and is usually referred to as congestion rent. 

Definition 3.8. The congestion rent on interconnector c = (r,s) e C at time t e T is given for n e 

l Axr andiel Cxr 



by: 



co c An, t) := (7i t - n c )t t 



Note that the congestion rent is indeed contained in the economic surplus in most market models 
run today. The surplus generated due to congestion rent is usually reinvested in order to reduce tariffs 
or build new lines j We will see later that the flows are not uniquely determined by the maximal 
economic surplus, but in general ambiguities only exist, when prices of adjacent areas coincide. In these 
situations, there is no price difference (tt s t — n T t ) and thus no congestion rent, and the congestion rent 
cannot be influenced by changing the flow t c t within the ambiguities. For more details on ambiguities 



of flows and how to choose unique flows refer to Section 5.3 



Definition 3.9. Let (tt, /3, ip, 5, t) be feasible for flMPECl ). The economic surplus of the market is the sum 
of the surpluses of all market participants and is computed as follows: 



■«Ih(5h) 



r rqiAOh) 

co(tz,I3,<p,5,t) = ^^ J] | (p h (q- a (u)) - % a>t ) 

aeAteT L heH a t , Aq h #(W0 



(3.15) 



beB a f£F a 

c=(r,s)eC t£T 



We will now show how the objective function can be simplified given condition (CI). The simplified 
expression of the economic surplus is equivalent to the objective of the welfare maximization problem 
in QNiu et al.| [20051 Section II.B]. 



Theorem 3.10. Let {n, /3, tp, 5, t) be within the given bounds and satisfy condition (C 1), then the eco- 
nomic surplus is given by 



f f 5h 

eo(/3, <p, 5) = J] J] S Aq M Ph^)dz +Y^Pb<lb,tPb+Y^ p f q f Vf.- 

aeAteT \heH at JO beB a feF a 



\ 



+ K 



= ^((Ph + Ap h ')Aq h 5 h --Ap h Aq h 5£)+^ i p b q bit p b + ^ i p f q f <p fit + K, 



h£H 



b£B 



f£F 
teT 



where K is a constant. 

Proof. By construction we have C = UaeA^o" = UaeA^a" • With this equation we get for all t S T 



2_j ( n s,t ~ n r,t) T c,t = 2_j n a,t I Zj Tc ' t ~~ Z-J Tc - C ) ' 
c=(r,s)eC aeA VeC" ceC+ J 



(3.16) 



Combining the above with constraint (CI) and substitution into fl3.15[ ) yields the desired result. □ 



Energinet.dk: 100% to reduce tariffs, Fingrid: 100% to build new lines, Statnett: 100% to reduction, Svenska Kraftnat: 
100% to build new lines, (cf. http://www.nordpoolspot . com/reports/Bottleneck- income/) 
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Observe that the objective function obtained in Theorem 3.10 is a quadratic, concave function and 
does not involve price variables n or flow variables t. While this is convenient for the actual com- 
putation, these variables still have to satisfy the stated conditions and so they implicitly influence the 
economic surplus. 

The problem formulation QMPECfl with the objective function as stated above is in principle solvable 
by convex Mixed-Integer Quadratic Programming solvers like IBM CPLEX if we linearize the price con- 
ditions by using the big-M method. Solving real world instances with this technique however is very 
hard, because each nonlinear price constraint must be modeled by using an auxiliary binary variable. 
For the considered instances we were not able to obtain feasible solutions or improve provided warm- 
start solutions within 30 minutes. This is clearly unsatisfactory as the actual market coupling auction 
needs to be cleared within 10 minutes. We therefore propose a heuristic and a specialized optimization 
algorithm for this problem. We did not yet try to solve the MPEC by using MINLP solvers like BARON 
or BONMIN, or general MIQP solvers like SCIP however we believe that due to the problem structure 
the presented problem specific approach will outperform general purpose solvers. 

4. Optimality Conditions 



In this section we define a relaxation of ( |MPEC| ) and analyze the optimality conditions of this prob- 



lem. Based on this analysis we derive a heuristic that performs very well in practice. We say that 
(/3*, i/?*) is a bid selection if §* e {0, 1} B and ip* e {0, l} FxT . Note that a bid selection does not neces- 
sarily satisfy the condition that none of the chosen bids should incur a loss. 

We consider the following relaxation of ( |MPEC| ) given a fixed bid selection fj8*, (/>*): 



f5 h 

max co(5, t) = Y Aq h p h (u)du + /(/3*, ip*) + K 

h£H JO 

s.t. [n] VatA,ttT q^ + J] Aq h 5 h + J] r c>t - J] T c>£ + k a>t (/3*, ip*) = 

heH a>t ceC+ ceC~ 

[v] Vfr e H 5 h -l<0 

[v] Vh e H ~5 h <0 (QPRelax) 



[p] 


Vc € C, t S T 


T c jt -^c,t <0 


[p] 


Vc e c, t e t 


tc,t - T Cjt _i - T c < 


M 


Vc<EC,t<ET 


-T C;t + T C;t <0 


[p] 


Vc e C, t e T 


T C ,t-l - ^c,t - ^ c < 



Note that ( |QPRelax] ) arises from ( |MPEC| ) by dropping the combinatorial structure, i.e., block bids 



and the flex bids. In this case we do not need to require the filling condition ( C 8| ) or the flow price 



condition ( |C 7\ explicitly as these naturally follows from the Karush-Kuhn-Tucker optimality conditions 
(cf. [Boyd and Vandenberghe, 2004, Chapter 5.5.3]). The next theorem shows how an optimal solution 
to ( |QPRelax| ) can be extended to a feasible solution of flMPEC[ ). 

Theorem4.1 (natural spot price characteristics). Let (/3*,i/?*) be a bid selection, f(P*,ip*), fc a t (/3*, c/?*) e 
K constants, and (5*, t*) an optimal solution to ( jQPRelax] ). Then there exist prices n, so that (5*, tt) sat- 
isfies the filling condition and (t*, n) satisfies the flow price condition. Furthermore, with 

f(fi*,<P*)= Xi Pb a b,tPt+ Xi Pf1fV*f,t aTld 

b£B,t£T feF.teT 

KtU 3 *' V*) = Yi qh >t ft + Yi q f V*f,t for at A, t£ T, 

b£B a feF a 
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the economic surplus in point (5*, t*) is given by the objective value co(5*, t*). 



Proof. Let (/3*, ip*) be a bid selection and let (5*, t*) be a maximum solution of ( QPRelax] ). As co is 
continuous differentiable and the constraints are linear, optimal solutions to ( QPRelax] ) have to satisfy 
the Karush-Kuhn-Tucker optimality conditions: There exist variables 71, v, v, ju, /i, p , p with 



Va<EA,t<E T,h<EH at 
yc = (r,s)<EC,t<E T 

Va<EA,t<E T,h<EH at 
\/c = (r,s)<EC,t<E T 



v,v,^,^,p,p >0 

-Aq h p ft (5*) + n a _ t Aq h + v h -v h = 

QX ct - M c f ) + iPc,t ~ Pc,t) + [(Pc,t+1 - Pc,t+l)] + n r,t ~ n s,t = 

(5 h -l)v h = 

(-5 h )v„ = 

(T C)t - T C)t ) JI c>t = 

( T c,t - T c,t-1 - T c ) p C]t = 



(t 



c,t-l 



'C,t 



T c )p t = 0, 



(4.1) 
(4.2) 
(4.3) 

(4.4) 
(4.5) 
(4.6) 
(4.7) 
(4.8) 

(4.9) 



where the term in brackets is only added if t + 1 e T. The conditions ( |4.1| ) to Q4.3D are the dual 
feasibility conditions and conditions ( |4.4| ) to ( |4.9| ) are the complementarity conditions. Observe that 



equations Q4.1D , Q4.3D , and ( |4.6D to ( |4.9D imply that (t*, tt) fulfills the flow price condition (C 7). The 



filling condition (C 8) can be derived from the equations ( |4.1D , Q4.2D , and ( |4.4| ) to ( |4.5| ). For this we 



consider three cases. Let a e A, t e T, and fr G H a t with Aq h ^ 0. In the first case, let 5^ = 1. Then 
v h = and v h > 0, and Q4.2J ) simplifies to — p^(5p + 7i a t < 0. For the second case suppose 5^ e (0, 1). 
Then v h = and v ft = 0, and ( |4.2| ) simplifies to — p^(5p + 7i a t = 0. For the last case, let 5* h = 0. Then 
v_ h > and v h = 0, and ( |4.2| ) simplifies to — p?,(5p + 7i a t > 0. Therefore (5*, 7i) satisfies the filling 
condition. 

For the last statement of the theorem, let/ and g a t be as above. As (5*,t*,/3*, (/?*) satisfies constraint 
(CI) we can apply Theorem 3.10 and obtain that the economic surplus is given by o>(<5*, t*). □ 



Observe that / and k a t in the proof above correspond to the contribution to the welfare (not neces- 



sarily positive) and the additional volume generated by the bid selection. We will later use Theorem 4.1 
to extend solutions to ( |QPRelax| ) together with a bid selection to feasible ones for flMPECl ). 



5. Uniqueness of a Solution 

In the following section we will discuss uniqueness aspects of the computed solution. More precisely, 
we will show that the economic surplus for a given bid selection is unique. The prices however do not 
necessarily be unique as we will see. 



Let7ieP AX ',5e[0,l] w ,T 



>CxT 



, and t c t < x ct < r ct . If n are prices so that (5, n) satisfies the 



filling condition and (t, n) satisfies the flow price condition, then we say that p are shadow prices for 
(5,t). 



5.1. Unique economic surplus for a given bid selection. Recall that the model ( |MPEC| ) requires the 
filling condition and the flow price condition to be satisfied. We can therefore only construct a solution 
to QMPEC1 ) from a solution to ( |QPRelax] ) if the solution satisfies the filling and the flow price condition. 
Fortunately, these are induced by the KKT optimality conditions of ( |QPRelax] ). Therefore every solution 
of ( [QPRelax"! ) for which shadow prices exist is optimal for ( [QPRelax] ). As there exists exactly one optimal 
objective value per bid selection, there exists exactly one economic surplus per bid selection. 
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Figure 5.1. Both flows generate the same economic surplus. 



5.2. Uniqueness of hourly bid execution. Let (/3*, if*) be a bid selection and let I c H be the indices 
of hourly bid segments with slope neither nor oo, i.e., 

I := {h e H | Aq h # and Ap h ^ 0} . 

An optimal solution to ( |QPRelax| ) is unique with respect to 5 7 , as co is strictly concave with respect to 
5 j (cf. |Dattorrol 02005] ). 

5.3. Choosing unique flows. The constraints in Section 3.1 and the optimality of a solution to ( |QPRelax| ) 
are not sufficient to imply a unique flow. This follows from the fact that even some of the vertical seg- 
ments (i.e., those with Ap h = 0) of hourly bids are not unique. To obtain a unique flow, we minimize 



the squared flow while fixing the economic surplus. In Figure 5.1 an example of a non-unique flow is 
depicted. The left part shows the solution with minimized squared flows. 

The following model selects a unique flow: Let (/3*, (/>*) be a bid selection, (5*, t*) be an optimal 
solution to ( |QPRelax] ), and let J c {h e H \ Aq h ^ and Ap h = 0} be the set of non-unique vertical 
segments. 



mm 



s.t. 



teT ceC 



(FixFlow) 



The economic surplus must not change (cf. Theorem 3.10) 

^p h Aq h 5 h = ^p h Aq h 5* h 



heJ 



heJ 



Clearing condition (cf. (C 11): 
VaeA,teT 

h£JnH at ceC+ ceC~ heJnH at ceC+ ceC~ 

Upper bounds, lower bounds, and ramp rate: 
Vh e J 0<5 h <l 

Vc e C, t e T Tct- T c,t - ^c,t 
Vc S C, t S T - t c < T ct - T,.^! < t c 

Note that the objective is strictly convex and the constraints are affine linear. So the obtained solution 
is unique with respect to the flow. The unique flow determines a unique net export per area and hour, 
therefore also the execution of vertical segments of the net curves is unique. 

It is easy to see, that ambiguities of the flow can arise in situations where the prices of adjacent areas 
coincide. As mentioned above there is no congestion rent in these situations. Our hypothesis is, that 
whenever there are flow ambiguities the prices of the involved areas coincide or can be chosen, so that 
they coincide. 
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5.4. Choosing unique prices. We will explain an easy way of how to obtain unique prices. Let (/3*, ip*) 
be a bid selection, (5*, t*) be an optimal solution to ( |QPRelax| ) with unique flows, and n be shadow 
prices. If in a specific area, in some time slot the net curve is horizontal (i.e., no volume change) at 
the clearing price, then the price can be chosen freely within the price indifferent price range without 
affecting the economic surplus (provided that the price does not violate any other constraint). We 
therefore propose, as above in the case of flows, to minimize the squared prices in order to determine 
unique shadow prices. 



2X 

aeAteT 



-2 
■a,t 



(QPPrice) 



s.t. filling condition with 5*, t*, /3*, ip* fixed 

flow price condition with 5*, t*, /3*, if" fixed 
executed block bids must not incur a loss: 

VaeA, be B a with /3* = 1 : 



^(Pb-^a,t)<lb,t >0 



ter 
executed flexible bids must not incur a loss: 

(Pf ~ n a,t) sgn(q / )>0 



Va eA,/ S F a , t e T with Vf t = 1 



Observe that the filling condition, the flow price condition, and the bid price constraints simplify to 
linear constraints, because the flow t*, the hourly bid execution 5*, and the bid selection (/?*, ip*) are 
fixed. 



6. The bid cut algorithms 

We will now present the proposed algorithm to solve the clearing problem. The basic idea is to relax 
all price constraints. For example we allow that block or flex bids incur a loss. The natural spot price 



characteristics (see Theorem 4.1 1 guarantee the existence of shadow prices for an optimum solution of 
the following relaxation: 



max 
s.t. 



co(j3,ip,5) 



CI): clearing condition 



C 2 1 : ramp rates 



C 3J ) : block bid links 

C 4 1 : flexible bid execution 

^£,e{0,l} forbeB 

V7 /it e{0,l} for/eF,ter 

<5 h e[0,l] for hen 

T c,t e [l c t^ctl forceC, teT 



(QPBidCut) 



Theorem 6.1. Let (/S*, ip*,5*, t*) be an optimum solution to ( QPBidCut ). Then there exist shadow 
prices nfor (5*, t*) and the objective value o>(/3*, ip* , 5*) coincides with the economic surplus in the point 
(£*,</>*, 5*, T*). 



Proof. Let (/S*, ip*,5*,z*) be optimal for ( IQPBidCut] ). Then (I3*,<p>*) represents a bid selection and 
(5*, t*) is optimal for ( [QPRelax] ). With Theorem 4.1 the existence of shadow prices is guaranteed and 
the objective value is the economic surplus. □ 
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Note that the theorem above only guarantees the existence of shadow prices, but does not guarantee 
that none of the executed block bids incurs a loss. 

Suppose now that we have an optimal solution to ( IQPBidCut] ) and the associated shadow prices. 
Then there might be block or flexible bids that incur a loss. Thus the solution is not valid for ( |MPEC| ). 
In this case an additional constraint, a so called bid cut, is added to the problem. This constraint 
cuts off the invalid bid selection and a class of similar bid selections. This process is repeated until 
a selection is found that is valid for flMPECl ). Shadow prices for a solution to ( IQPBidCut] ) can be 
computed with the model ( |QPPrice| ), whereas the bid price constraints need to be relaxed to make 
( JQPPrice] ) always feasible. This can be done for example by introducing positive continuous variables 
A b for every executed combinatorial bid b and adding X b to the left side of the according bid price 
constraint. Then we add all these variables to the objective function, so that they get minimized. Now 
the variable X h models the loss of the executed combinatorial bid b. 

We will now introduce two families of cuts that can be used. The first one works very well in practice 
but is slightly inexact as it might converge to a slightly suboptimal solution. The second family of cuts 
is exact but more cuts have to be added. We will later report timings for both families in Section [7J 



Definition 6.2. Let (/3*,c/?*,<5*,t*) be a solution of ( IQPBidCut] ) and n* be shadow prices for (<5*, t*). 
The set of bids, that incur a loss at prices n*, is given by 

B loss = {b^B\l3* b = land JJp b - n* at )q b>t < 0} 
and F loss = {f e F | ex. t e T with <p* >t = 1 and {p f - tt* ^sgn^) < 0}. 



were B loss is the set of executed block bids that incur a loss at prices n* and F; 05S is the set of executed flex 
bids that incur a loss at prices ti*. The linear constraint that prohibits the execution of at least one bid of 
the set B loss U F loss is given by 



Cut(B; 0SS ,F; 0SS ) : 



2> 



+ 



b£B, r 



f&h 



,teT: V * fjt * 



<Pf,i 



< \Bi n ,A + I F; n r J — 1 



lossl 



We refer to this class of constraints as Bid Cuts (cf [Mutter 2009 Definition 5.2]). 



As mentioned above, for this family of cuts, the obtained solution might be slightly suboptimal, i.e., 
these cuts provide us with a heuristic. Nonetheless, this heuristic can be used iteratively as shown 



in Algorithm 6.1 or it can be incorporated into a branch-and-cut framework. The latter is the more 
effective one as it provides us with an exact solution algorithm. 

If the Bid Cut is replaced with a less aggressive one that only removes exactly one invalid bid selec- 
tion, then the resulting algorithm will converge to the globally optimal solution. This Branch-And-Cut 
Decomposition, works, for example, with the following family of cuts that separate exactly one invalid 
bid selection. 



Definition 6.3. Let (/3*,(/>*,5*,t*) be a solution of ( IQPBidCut] ). Then following inequality cuts off the 
bid selection (/J*, (/>*): 

Cut(0>*): X! h+ S VU + £ d-foH S (1 -¥>/,*)>! 

beB:Pl=0 f£F,t£T-.Lp* =0 beB:/3*=l f£F,t£T:ip* =1 



We refer to this class of constraints as no-good cuts (cf. ftMuller 2009 Theorem 5.4], Chu and Xia 



Note that the combinatorial bid curtailment (Remark[T]) must be ensured in a similar way. We cannot 
add it directly to the ( QPBidCut] ) -model, as it would affect the properties of an optimal solution. 

20 



Algorithm 6.1 Iterative Bid Cut Heuristic 

Input: An instance of the model ( |QPBidCu"tj ) 
Output: (re*, 13*, (/>*, 5*, t*) feasible for QMPECD 



done <— false 
while ->done do 



(P*,(p*,5*,T*) «- Solve flQPBidCutp . 



(t*) <— Solve (FixFlow) with input values (/?*, ip*, 5*). 
in*,B[ 0SS ,Fi 0SS ) «- Solve QQPPriceD with input values (/T, </>*, 5*, t*). 



if l B los S U^ ss l = 0,then 

done <— true 
else 

Add Cut(BJ 0SS ,F^ ss ) to the model ( |QPBidCut| ). 
end if 

i«-i + l 
end while 

Return (7r*,/3*,</>*, 5*, t*). 



7. Results 

We will now present some computational results for the algorithms presented in Section [6j We 
compare the Bid Cut Heuristic and Branch-and-Cut Decomposition to the commercial algorithm EMCC 
Optimizer of Deutsche Borse Systems which is currently used by European Market Coupling Company 
(EMCC) to determine flows between Denmark (Nord Pool Spot AS) and Germany (EPEX Spot). 

For our computational tests we used 79 realistic instances, that contain 10 European market areas, 
about 600 combinatorial bids, and about 31.700 defining net curve nodes. All tests were performed 
on the same hardwardj using IBM CPLEX 12.1 as mixed integer linear/quadratic programming solver. 
In our tests, we did not incorporate the ( |FixFlow| ) model, as its sole purpose is the redistribution of 
flow without affecting the overall economic surplus. Moreover, our tests indicated that fixing the flow 
impacts the prices only in very pathological cases that can usually be disregarded. We will first give a 
brief introduction to the current version of the EMCC Optimizer and provide some statistics. Then we 
discuss the two proposed algorithms from above. 

EMCC Optimizer. The version of the EMCC Optimizer that we used in our comparison is version 2.1.2 
(in use since 09.11.2009). This algorithm is based on computing a start solution with a linear program 



(cf. Krion [2008]) and then checking the feasibility of the bid selection. If the bid selection is not 



feasible, executed bids are excluded successively until the solution is feasible. Then the solution will be 



improved by trying to include not executed bids. For more details see [EMCC GmbH Chapter 4.3.2]. In 
some cases this search phase is very time consuming and cannot be completed within 10 minuteqj In 
this case, the last feasible solution will be returned; in fact the start solution for the last phase is already 



7 Intel Xeon Core E5440 with 8 GB memory 

The limit of 10 minutes reflected the operational requirements when the tests of this algorithm were performed (Septem- 
ber 2009). For production it has been changed to 15 minutes. 
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Algorithm 


Optimal Bid 
Selection 


Best known 
Bid Selection 


PRBs 


Relative 
Gap 


Computing Time 


Stopped at 
max Time 


Average (s) 


Min (ms) 


Max (s) 


EMCC Optim. 2.1.2 


> 15% 


24% 


10.76 


2.577E-6 


361.4 


46000 


600.0 


22% 


Relax And Fix 1 


> 35% 


71% 


5.04 


2.01 3E-6 


0.7* 


193* 


1.6* 


0% 


Bid Cut Heuristic 


> 38% 


96% 


4.65 


1.926E-6 


4.1 


493 


62.8 


0% 


B&C Decomposition 


> 38% 


1 00% 


4.62 


1.924E-6 


(9.3) 390.8 


1114 


600.0 


65% 



1 The test of this algorithm was performed on a 1.8 times faster machine, a Core 17 920. 

Table 1. Numerical results for 79 realistic instances 



feasible however might be suboptimal. In our analysis, for 22% of the instances this last heuristic phase 
could not be completed on time. The average time for computing clearing prices was 6 minutes and 
in at least 15% of the cases the bid selection was optimal. In 76% of the cases the selections could be 
slightly improved by our algorithms and in the remaining cases we could neither prove nor disprove 
the optimality of the selection. The relative distance to the best upper bound, i.e., the relative gap, 
averages to 2.577xl0~ 6 ; a selection is considered optimal here if the gap is below lxlO -12 which 
roughly corresponds to one euro cent. On average there are 10.76 rejected combinatorial bids, so 



called paradoxically rejected bids (PRB) (cf. Meeus et al. [ 2009 ]), that could potentially have generated 
a net profit. However, in most cases these bids cannot be included without decreasing the overall 
economic surplus or rejecting other bids that are in the money. 

Bid Cut Heuristic. The tests are based on a branch-and-cut version of the Bid Cut Heuristic. In at least 
38% of the tests the bid selection is optimal and the gap averages 1.926x 10~ 6 . Only in 4% of the cases, 
the solution could be improved and in the average there are only 4.65 PRBs leading to a potentially 
higher overall acceptance of clearing prices/decisions. The most important property of the algorithm 
however is its running time which averages to only 4.1 seconds. The maximum computing time was 1.1 
minute and so clearing within the 10 minutes is easy possible. Moreover, due to being able to complete 
the process before the time limit is reached the results are reproducible on other machine^ 

B&C Decomposition. The Branch-and-Cut Decomposition is slower than the Bid Cut Heuristic. How- 
ever using it, we can prove optimality of a solution. It uses the solutions of the Bid Cut Heuristic as 
warm start solutions. In 38% of the cases the algorithm finishes in about 9 seconds. This implies that 
the found bid selection is optimal in at least 38% of the test. In the other cases the solver was stopped 
when reaching the time limit without proving optimality. During the computation 4% of the bid selec- 
tions found by the Bid Cut Heuristic could be improved. This reduced the relative gap to 1.924xl0 -6 
and the number of PRBs was reduced to 4.62. The improvement which is not significant indicates that 
the solutions obtained via Bid Cut Heuristic are very good. 

In Table [l] we summarize the results. More details about each instance can be found in Tables |2] and 
[3} We can see, that the relative gaps of all algorithms have the same order of magnitude. The best 
solutions are found by the Branch-and-Bound Decomposition, but the relative gap cannot be reduced 
significantly. The fastest algorithm is the Bid Cut Heuristic, that finds very good bid selections in about 
4 seconds. 



8. Concluding remarks 

Our proposed approach works very well in practice and we were able to derive nice properties from 
the formulation as an optimization problem. A possible route for improvement is to replace the bid cut 



by an infeasibility cut of the generalized Benders decomposition (see Geoffrion [ 1972 ], Chu and Xia 



Reproducible prices are crucial for market transparency and confidence 
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Table 2. Detailed results, test cases 1-40. Computing times are measured in milliseconds. 



# 


Comb. 


Net Curve 


Solver 2.1.2 




Bid Cut Heuristic 




BnC Decomposition 






Bids 


Segments 


Comp. Time 


Rel. Gap 


PRBs 


Comp. Time 


Rel. Gap 


PRBs 


Comp. Time 


Rel. Gap 


PRBs 


1 


742 


34536 


299000 





4 


2632 





4 


6621 





4 


2 


727 


34726 


* 821000 


4.00E-10 


8 


3101 





7 


600565 





7 


3 


792 


35227 


208000 


3.19E-09 


5 


2017 


3.19E-09 


5 


600562 


3.19E-9 


5 


4 


750 


34733 


476000 


4.77E-08 


10 


2621 


4.61 E-08 


7 


600334 


4.61 E-8 


7 


5 


685 


32613 


190000 


4.64E-08 


10 


1881 





6 


4929 





6 


6 


630 


30065 


350000 





3 


1086 





3 


3992 





3 


7 


775 


34516 


* 883000 


1 .34E-08 


8 


2615 


1.31 E-08 


7 


600251 


1.31 E-8 


7 


8 


772 


33545 


243000 


2.22E-07 


8 


2260 


1.70E-07 


6 


600339 


1.70E-7 


6 


9 


542 


34876 


578000 


2.30E-07 


5 


2528 





1 


12653 





1 


10 


595 


31226 


576000 


1 .99E-08 


6 


2199 





5 


53924 





5 


11 


585 


31888 


116000 


2.51 E-06 


17 


10102 


9.49E-07 


5 


600247 


8.85E-7 


5 


12 


598 


33262 


396000 


4.51 E-07 


8 


1968 





3 


9831 





3 


13 


435 


29590 


112000 


2.74E-06 


9 


9970 


1.14E-06 


7 


600213 


1.14E-6 


7 


14 


558 


31959 


460000 


3.63E-06 


18 


4531 


1.84E-06 


5 


600232 


1 .84E-6 


5 


15 


593 


35384 


399000 





2 


2143 





2 


7845 





2 


16 


583 


36302 


248000 


3.29E-06 


9 


2868 


2.83E-06 


2 


600339 


2.83E-6 


2 


17 


585 


37046 


602000 


3.78E-07 


11 


5676 


1.88E-07 


5 


600272 


1 .88E-7 


5 


18 


155 


18195 


111000 





1 


493 





1 


1114 





1 


19 


149 


18287 


61000 





1 


543 





1 


1135 





1 


20 


453 


30471 


141000 


2.29E-07 


9 


1235 


2.20E-07 


2 


600178 


2.20E-7 


2 


21 


201 


30027 


123000 


3.21 E-06 


12 


1871 





3 


6824 





3 


22 


675 


30027 


149000 


1 .73E-06 


12 


1824 


1.50E-06 


5 


600225 


1.50E-6 


5 


23 


701 


31559 


135000 


5.85E-06 


27 


34961 


5.92E-07 


12 


600210 


5.92E-7 


12 


24 


728 


29902 


440000 


3.43E-06 


4 


1334 


3.43E-06 


4 


600183 


3.43E-6 


4 


25 


651 


30196 


268000 


3.24E-06 


11 


2366 


3.11 E-06 


7 


600216 


3.11E-6 


7 


26 


642 


28589 


162000 


5.18E-06 


4 


2090 


5.16E-06 


3 


600237 


5.16E-6 


3 


27 


622 


31692 


342000 


3.40E-06 


21 


2720 


3.25E-06 


7 


600177 


3.25E-6 


7 


28 


604 


31535 


268000 


1 .08E-07 


4 


1551 


2.59E-08 


3 


600204 


2.59E-8 


3 


29 


707 


33882 


611000 


5.34E-06 


7 


3430 


5.33E-06 


6 


600201 


5.33E-6 


6 


30 


665 


33821 


249000 


1 .04E-08 


2 


1656 


1.04E-08 


2 


600285 


1 .04E-8 


2 


31 


684 


33975 


609000 


2.39E-05 


68 


62814 


5.04E-06 


12 


600231 


5.04E-6 


12 


32 


599 


27453 


181000 


8.31 E-06 


6 


2597 


8.32E-06 


6 


600243 


8.31 E-6 


6 


33 


420 


27247 


176000 


9.49E-06 


3 


1227 


9.49E-06 


3 


600222 


9.49E-6 


3 


34 


163 


24801 


46000 








802 








2620 








35 


142 


26142 


99000 


2.40E-08 


1 


1020 





1 


3090 





1 


36 


618 


31584 


362000 


1 .57E-07 


12 


2189 


2.23E-07 


7 


600226 


1 .23E-7 


7 


37 


249 


27231 


94000 





1 


1003 





1 


4403 





1 


38 


609 


29589 


152000 


3.78E-07 


27 


2279 


4.41 E-09 


11 


600224 


4.41 E-9 


11 


39 


297 


30596 


102000 


7.50E-07 


7 


1342 





1 


5729 





1 


40 


318 


30068 


247000 





2 


1681 





2 


6123 





2 



) Computing time was extended in this test, so that the inclusion phase could be finished. 
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Table 3. Detailed results, test cases 41-79. Computing times are measured in milliseconds. 



# 


Comb. 


Net Curve 


Solver 2.1.2 




Bid Cut Heuristic 




BnC Decomposition 






Bids 


Segments 


Comp. Time 


Rel. Gap 


PRBs 


Comp. Time 


Rel. Gap 


PRBs 


Comp. Time 


Rel. Gap 


PRBs 


41 


253 


27288 


158000 


6.43E-08 


5 


1062 





2 


4626 





2 


42 


604 


29969 


358000 


7.12E-05 


2 


1459 


7.12E-05 


2 


600261 


7.12E-5 


2 


43 


644 


32188 


361000 


1.96E-07 


9 


1658 


2.89E-08 


4 


600246 


2.89E-8 


4 


44 


683 


32477 


232000 





2 


1387 





2 


4125 





2 


45 


695 


33131 


608000 


1.64E-06 


7 


6630 


1.62E-06 


4 


600293 


1.62E-6 


4 


46 


576 


28686 


142000 


1 .07E-06 


22 


6151 


3.83E-07 


7 


600215 


3.83E-7 


7 


47 


476 


26189 


181000 


1.27E-09 


2 


752 





1 


1862 





1 


48 


655 


29370 


384000 


3.90E-08 


10 


1666 





8 


16594 





8 


49 


673 


32226 


396000 


1 .37E-05 


12 


3937 


1 .34E-05 


7 


600198 


1.34E-5 


7 


50 


667 


33513 


349000 


6.50E-07 


11 


1817 





1 


8442 





1 


51 


738 


32883 


621000 


2.36E-07 


12 


2499 


3.08E-08 


5 


600220 


3.08E-8 


5 


52 


628 


29959 


243000 


4.84E-07 


26 


1508 


2.96E-07 


3 


600195 


2.96E-7 


3 


53 


439 


29670 


324000 





6 


1607 





6 


3684 





6 


54 


705 


33437 


604000 


1.66E-06 


18 


2641 


2.72E-09 


5 


600240 


2.72E-9 


5 


55 


691 


33609 


623000 


3.53E-07 


10 


2148 





4 


6236 





4 


56 


696 


34061 


494000 


6.18E-08 


6 


2210 


2.70E-08 


4 


600250 


2.70E-8 


4 


57 


692 


33364 


158000 


9.78E-07 


14 


3282 


7.97E-07 


5 


600246 


7.97E-7 


5 


58 


644 


31157 


254000 


9.13E-07 


15 


3023 


7.30E-07 


6 


600230 


7.30E-7 


6 


59 


562 


30123 


303000 


1.05E-08 


5 


1388 





4 


5158 





4 


60 


603 


31833 


567000 


2.04E-09 


7 


2545 


2.04E-09 


7 


600225 


2.05E-9 


7 


61 


660 


34087 


620000 


6.46E-07 


18 


2259 


1.89E-07 


4 


600293 


1.89E-7 


4 


62 


682 


35228 


611000 


5.73E-08 


14 


2229 


3.25E-09 


5 


600281 


3.25E-9 


5 


63 


675 


35128 


345000 





2 


1561 





2 


5709 





2 


64 


713 


34921 


362000 


1.62E-07 


8 


2180 


1.30E-07 


5 


600326 


1.30E-7 


5 


65 


627 


32606 


203000 


1.78E-06 


21 


1607 


2.79E-08 


5 


600205 


2.79E-8 


5 


66 


680 


35368 


389000 





5 


2090 





4 


7680 





4 


67 


674 


34825 


298000 


6.11E-06 


12 


7555 


5.75E-06 


8 


601153 


5.75E-6 


8 


68 


711 


33534 


610000 


4.63E-07 


24 


1981 





5 


5924 





5 


69 


730 


32839 


626000 


1.12E-06 


20 


5031 


7.99E-07 


9 


600253 


7.99E-7 


9 


70 


697 


31728 


341000 


6.98E-07 


23 


1954 


4.81 E-08 


4 


600186 


4.81 E-8 


4 


71 


737 


33662 


319000 


2.36E-07 


22 


2816 


8.64E-09 


10 


600230 


8.64E-9 


10 


72 


708 


34162 


497000 


4.22E-07 


9 


1850 


3.43E-07 


6 


600247 


3.43E-7 


6 


73 


765 


34315 


403000 


7.74E-07 


7 


4061 


4.81 E-07 


7 


600262 


5.08E-8 


5 


74 


758 


36456 


* 745000 


6.83E-07 


18 


3424 


5.15E-07 


7 


600287 


5.15E-7 


7 


75 


542 


34437 


631000 


1.49E-07 


4 


2029 





1 


7967 





1 


76 


595 


30812 


610000 


3.31 E-06 


18 


2509 





3 


600227 





3 


77 


585 


31593 


410000 


1 .57E-06 


13 


13883 


1.30E-06 


6 


600248 


1.30E-6 


6 


78 


598 


32966 


479000 


4.80E-07 


12 


2093 





3 


51173 





3 


79 


558 


31525 


605000 


3.87E-06 


16 


21331 


1 .52E-06 


6 


600212 


1 .52E-6 


6 


avg. 


595 


31692 


361380 


2.58E-06 


10.76 


4114 


1.93 E-06 


4.65 


390808 


1.92E-6 


4.62 



) Computing time was extended in this test, so that the inclusion phase could be finished. 
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[2004]). The main difficulty here is to find a non-trivial cut that is strong enough to work as fast as 



the bid cut. Otherwise the performance of the algorithm for day-ahead market clearing might not be 



sufficient enough. For example the simple no-good cut (cf. Definition 6.3 1 used for the Branch-and-Cut 
Decomposition is a valid Benders infeasibility cut, but the resulting algorithm is to slow. 

A nice additional consequence of the presented model are the transparent pricing rules. In fact, for 
a given bid selection optimality can be checked easily: A market participant could collect all relevant 
information to check the flow price condition. It is clear, that a solution will also satisfy the filling 



condition, so by applying theorem 4.1 the participant knows that the economic surplus is maximal for 
the actual bid selection. In this case the market is in an equilibrium. This transparency is crucial for 
market confidence and liquidity. 

We want to mention, that the flow price condition is subject to change in the future. This is due to 
the fact that we made the assumption that no energy is lost during transmission. Unfortunately this 
assumptions is not realistic for some interconnectors. The loss of electricity, the so called deadband, 
has to be incorporated in the next model. This modification will alter the optimality conditions of the 
connectors that have a deadband. 
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order book 


participant 


a 


b 


c 


d 






price limit (€/MW) 
quantity (MW) 


1 
-1 


2 
1 


3 

-2 


4 

2 




solutions 


participant 


a 


b 


c 


d 


objective 


A discriminatory 


executed 

variable payment (€/MW) 

fixed payment (€) 

total payment (€) 

avg. payment (€/MW) 

gain (€) 


yes 
3 
2 
-1 
1 



yes 
3 
-1 
2 
2 



yes 
3 

-6 
3 



yes 
3 
2 
8 
4 



3 


B nondiscriminatory 


executed 

variable payment (€/MW) 
total payment (€) 
gain (€) 


3 




3 




yes 

3 

-6 




yes 
3 
6 
2 


2 


C nondiscriminatory 


executed 

variable payment (€/MW) 
total payment (€) 
gain (€) 


yes 

1 

-1 




yes 
1 
1 
1 


1 




1 




1 


D nondiscriminatory 


executed 

variable payment (€/MW) 
total payment (€) 
gain (€) 


























Table 4. combinatorial auction example 



Appendix A. discriminatory prices versus nondiscriminatory prices 

We want to illustrate the difference between our nondiscriminatory prices and the discriminatory 
prices proposed in O'Neill et al. [ 2005 ] by analyzing the small example given in Table [4| There are 
four indivisible orders in the market. Two supply orders (negative quantity) and two demand orders 
(positive quantity). Participant a is willing to sell one unit if he receives at least one euro, participant 
b is willing to buy one unit if the price is lower than two euros, and so forth. Solution A is the one that 



is proposed by O'Neill et al. [ 2005 ] whereas our market model would choose solution B. 



In solution A there exists no linear price. This is the reason why it is suggested to divide the price into 
two parts. A variable payment depending on the quantity and a fixed payment. A negative payment 
indicates that the participant receives a payment. As every participant pays or receives an individual 
fixed payment we see that this is a discriminatory pricing scheme. Furthermore we observe, that no 
participant is realizing any gain in this solution. If we look at the total payments that the participants 
pay to the auctioneer or receive from the auctioneer we observe that the auctioneer keeps the entire 
economic surplus that was generated by the indivisible orders. 

In solution B we simply do not execute the orders of participants a and b. This allows us to find 
linear prices, that are consistent with the dual prices that we know from the associated continuous 
optimization problems. The linear prices and the absence of a bid-ask spreads guarantee that the auc- 
tioneer cannot pocket any surplus. Indeed, the entire economic surplus is received by the participants. 
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Figure B.l. Active ramping condition in forward direction, ramp rate = 500 MW 

Appendix B. Simplifications of the flow price condition 
Proof of eq. ( |3.7| ) . Suppose that there is no ramping on any interconnector. Then we have t c = oo for 



all c e C, so that with fl3.5] ) and Q3T6J) it follows, that p ct = p ct = for all c e C and t e T. Equation 
( |3.2| ) before simplifies to 



(Met - /O + 



71,. 



»V = o. 



(B.l) 



We also know, that t c t < t c t implies /i c t = and that t c t > t c £ implies // =0. Together with (JbTTJ) 
we therefore conclude 



Vc = (r,5)eC, t S T 



<0 



7T C f — 71 



r.i ~ (Mc.t M c t ) | > Q 



'lc,t <T Ct 



Using this equation we can interpret the term Jl c t — /i as the shadow price of available transmission 
capacity: an additional unit of transmission capacity increases the economic surplus by \JI C t — \i 
euros (on a sufficiently small interval). 



□ 



Proof of eq. ( |3.8| ). Suppose that there is ramping on the interconnector c e Cand that the ramping 
condition (C 2) is only binding in the two consecutive time slots t and t + 1. Suppose further that 



there is still available transfer capacity in that time window. The available transfer capacity implies that 
Ic,t < T c,t < ^c,t and T c>t+1 < T C)t+1 < T ct+1 , so that ju Cjt = M c t = i" Cj t+i = M c t+1 = follows using 
( |3.3| ) and ( |3.4[ ). We assume that the ramping condition is active only between hour t and t + 1, and 



that it is active in forward direction (see Figure B.l I. This yields that — t c < (t c t+1 — t c t ) = t c and for 
all u e T with u ^ t we have — t c < (t c u — t cu _ 1 ) < t c . Thus p cu and p cu vanish, except for p c t+1 . 
From equation ( |3.2[ ) and ( |3.1D we get 



Kr,t- n s,t=Pc,t+l ^0 and 
^r,t+l _ rt s,t+l = - Pc,t+1 < 0. 

These statements together provide — (7t r t — n st ) = — p c t+i = ( n r,t+i ~ t s t+i)j ^r t > ^s.t' an< i ^rt+i < 
7T S t+1 . We can now interpret the variables p as the shadow prices for ramping, in short the ramping 
prices. Here the economic surplus could by increased by p ct+ \ euros if the ramp rate would be increased 
by one unit (again on a sufficiently small interval). □ 



Proof of eq. fl3T9] ) and ( [3.10D . See also [ ]Muiler[ [2009] Theorem 4.4]. Let the flow-price pair (t, tt) 
satisfy the filling condition, let c = (r,s) G C, teT 
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and 7i rt < 7T S t . 



Assume that r ct < T ct and 



^c.t - T c,t-i < T c [. and T c,t - T c,t+i < T cL then we have ju ct = 0, p ct = [, and p ct+1 = 0]. 
Consequently we obtain with fl3.2[ ) and ( |3.1| ) 



( Mc,t -/*) + ( Pc,t - Pc,t ) + [Pc,t+1 - Pc.t+ll + ^r.t - ^s,t = ° 



>o _u >o 



>0 <0 



which is a contradiction. Implication ( |3.10P can be proven in a similar fashion. □ 



Appendix C. MIP formulation of the filling condition 

The filling condition can be easily formulated using binary auxiliary variables. For every segment 
h e H we introduce the binary variable y h e {0,1}. Without loss of generality we assume that the 
segments are sorted at first by areas and hours, then by ascending prices, and finally by descending 
quantities. The filling condition can now be formulated like this: 

Va e A, t e T, h e H at \ maxH at 5 h <j h < 5 h+1 . 

Note that the presented algorithms do not rely on the inclusion of this condition. In fact, including it 
does not impede the solving process. The slowdown arises due to the inclusion of the filling condition 
combined with the price conditions of combinatorial bids. This is due to having also to formulate the 
flow price condition by using a complicated MIP formulation and this condition cannot be solved fast 
enough by standard solvers. 

Appendix D. Presolving 

Before starting the optimization process, we can compute upper and lower bounds for the prices. For 
every area and hour we can determine the intersection of the hourly net curve with the price axis. This 
price would be the clearing price if no block bid would be executed and no power would be imported 
or exported. Now we assume that we have to import as much additional power as possible and that 
all supply block bids are executed. The additional power must be bought by hourly bids, thus the price 
decreases along the hourly net curve until all additional power is consumed by hourly bids. This price 
represents a lower bound for the specific area and hour. Upper bounds can be determined similarly. 



Formula ( |D.1| ) is directly derived from the clearing constraint (CI) and shows the upper and lower 



quantity bounds of the hourly net curves. The price bounds can be derived from these quantity bounds. 
q + Y Aq h 5 /- "^^ min{ °' qb - t} ~ ^ /efa min{0, q/ } +I] ceC - T c>t - 2 ceC T^ t ^ ^ 

With the help of these bounds the execution state of the blocks and hourly bids that always realize a 
gain or incur a loss is given. The bounds can also be used to cut off infeasible classes of bid selections 
in the branch-and-cut tree. We implemented theses techniques, but unfortunately the solving process 
could not be enhanced significantly. 
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